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Abstract 

A general structure is developed from which a system of integrable partial difference 
equations is derived generalising the lattice KdV equation. The construction is based 
on an infinite matrix scheme with as key ingredient a (formal) elliptic Cauchy kernel. 
The consistency and integrability of the lattice system is discussed as well as special 
solutions and associated continuum equations. 

1 Introduction 

In recent years the integrability of discrete equations, i.e. ordinary or partial difference 
equations as well as analytic difference equations, has become an issue of considerable 
attention, (cf. the biannual sequel of SIDE meetings on Symmetries and Integrability 
of Difference Equations, cf. http://www.maths.leeds.ac.uk/~side). A particularly 
interesting development has been on the one hand the search for precise definitions of 
integrability in the discrete regime and the subsequent development of integrability detec- 
tors, on the other hand progress in the classification of various types of integrable discrete 
systems, such as the so-called discrete Painleve equations. 

One line of research since the mid-seventees with the pioneering works of Ablowitz and 
Ladik, Q and of Hirota, |J, has been the study of integrable lattice equations, which 
are partial difference equations living on two-or higher dimensional space-time lattices. 
A systematic approach for the deriviation of such equations was performed in a series of 
papers, cf. e.g. fll7| , p3| , [2^| from the point of view of formal integral equations. The 
key ingredients in this approach, loosely coined direct linearisation, are two-fold: linear 
discrete or continuous dynamics residing in plane wave factors (the object p(k) introduced 
below), and a formal Cauchy kernel (0(«i,K2) m the notation of the present paper) by 
which the connection with nonlinear objects takes place. The implementation of the 
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method is in terms of an infnite matrix structure, which, albeit stricly formal, has proven 
to be powerful tool in setting up the basic relations from which all ingredients (Lax pairs, 
Backlund and Miura transformations, hierarchies of commuting flows, etc.) can be derived 
without having to rely on further Ansdtze. 

The resulting discrete systems on the 2D lattice are particularly important for devel- 
oping an understanding of the notions of discrete integrability, in particular since other 
reduced systems (such as discrete Painleve equations, integrable finite-dimensional map- 
pings and discrete-time many-body systems) can be viewed as special solutions. Further- 
more, there exist compelling connections with difference geometry, and the combinatorics 
on graphs. 

The simplest example is possibly the lattice (potential) KdV equation 



a + b + u — u) (a 



b + u 



(1.1) 



which, together with its companion equations and reductions, was studied at length in 
numerous papers, cf. e.g. |L7], 22, |2l], [7|, [Eg] . 

The notation used to describe this lattice system is as follows: Denoting the dependent 
variable on the two-dimensional lattice by by u = u n ^ m (n, m G Z) the equation is given 
in terms of the variables u, u, u and u around an elementary plaquette as indicated in 
Figure 1. 

U n U 



u a u 

Figure 1: Elementary quadrilateral on which the lattice equation is defined. 



The integrability of such lattice equations can be understood in a rather simple way: 
it seems entirely to reside in a simple but deep combinatorial property, first described 
in the paper [19|, namely that two-dimensional lattice equations in fact form parameter- 
families of compatible equations which can be consistently embedded in a multidimensional 
lattice, on each two-dimensional sublattice of which a copy of the lattice equation can be 
defined. As was shown in AO], 13], cf. also M, this property is powerful enough to derive 



subsequently Lax pairs for the lattice equations, which can then be used to study the 
analytic properties of solutions. More recently in (3| the property was used to arrive at a 
full classification of lattice equations integrable in this sense under certain some additional 
simplifying assumptions. In particular, the most parameter-rich equation emerging from 
this classification was a lattice system already found by V. Adler in Q , for which the Lax 
pair was derived in fl3[ |. 

Although these developments form a considerable step forward, relatively little is known 
of the underlying algebraic and analytic structures of Adler's lattice equations (this being 
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the subject of a paper []15[ in preparation). In this note we present another integrable 
lattice system which is naturally associated with an elliptic curve, and which is given by 
the following coupled system of equations: 



in which 



a + b + u — u) (a — b + u — u) = a 2 — b 2 + f (s 



s) s — s 



s — s (w — w) 



[(a + u)s — (6 + u)s] s 



[a — ujs 



(b — u)s 



(s — s) ( w — w j = (a — u)s + (b + u)s s— (a + u)s + (b — u)s 



a + u 
b + u- 



w 
s 
w 
s 



~ W \ 2 T^/ ~N 

a — u-\ J = a — F{ss) 



b-u + 



w 



P(ss) 



(1.2a) 
(1.2b) 
(1.2c) 

(1.2d) 
(1.2e) 



P(x) = - + 3e + fx , 

x 

with e and / moduli of an elliptic curve y 2 = P(x) . We will show that this system 
( |1.2j ) is a natural extension of the lattice KdV equation (1.1), the latter being obviously 
recovered when / = 0, i.e. when the elliptic curve degenerates. 



The system (1.2) is admittedly rather complicated, and is outside the class of systems 
considered in the paper Q. Rather than starting with an ad-hoc Ansatz we want to 
derive the system in a systematic way from a structure that not only provides us with 
the system but also allows us to construct its Lax pair and its linearisation. In order to 
achieve this we start from an infinite matrix system with a quasigradation of elliptic type 
similar to the ones that occur in the Krichever-Novikov algebras. The nonlinear structure 
arises from objects in this algebra defined through an elliptic Cauchy kernel. As a result 
we find a system of equations which forms a two-parameter "elliptic" deformation of the 



lattice systems of KdV type investigated in the past, cf. e.g. [17, 23, 14|. The extended 
lattice systems then will serve as a starting point for the study of a number of new discrete 
systems, notably the ones that arise from finite-dimensional or similarity reductions. Also 
their continuum counterparts form an interesting new class of integrable systems that 
merit further investigations. 



2 Infinite Matrix Structure 
2.1 Elliptic Matrices 

The starting point for the construction of the elliptic discrete model is an algebra A = 
Matg(C) of centred infinite "elliptic" matrices. We prefer to describe this algebra in simple 
terms without going into details of the construction, as we are only interested in some basic 
relations. 

We mean by A the associative algebra (with unit 1) with the usual matrix multiplication 
relative to a designated central entry, i.e. elements A,B G A are infinite matrices A = 
{Aij), B = (Bij), € Z), having ^4q,0! respectively i?o,o as central entry with respect 
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to which the matrix multiplication is defined, i.e. through the usual formula 
(A ■ B)ij = A iik B k j 

Since we will only be interested in the formal structure leading eventually to nonlinear 
discrete equations, we will not be distracted by questions related to the convergence of 
infinite sums (in which case one has to impose some additional restrictions on the matrices, 
which, however, will not be relevant for the present discussion). The main property of 
this elliptic matrix algebra is that it is not graded in the usual way (i.e. through the 
organisation of rows and columns according to integer steps), but quasigraded by means 
of two types of index raising operators A (of order 1) and P (of order 2), subject to the 
algebraic relation 

A 2 = p + 3el + IT' 1 , P A = A P, (2.1) 



where e, / £ C are moduli of an elliptic curve given by the relation fl2.l|) . This means that 
the indices of the matrices are defined relative to the action of these operators, namely as 
follows: for each element A = (Ay) £ A we identify the indices as follows: 

A 2it2j = (P* ■ A ■ *pi) 00 
^ + i, 2j = (P l -A.A.^') 00 

A 2i>2j+1 = {P i -A- t A- t P% (2.2) 
A 2i+ i, 2j+ i = (-P* A - A - ' A • 'P^ 00 

We need three ingredients to formulate the scheme: 

• A parameter-family of elements C £ A which encodes the dynamics of the system 
via the linear relations 

C ■ (a - *A) = (a + A) • C , (2.3a) 
C-(6-*A) = (6 + A)-C. (2.3b) 



In addition we can impose differential relations with respect to a hierarchy of con- 
tinuous flow- variables xf 

-^-C = A? • C - C(- *AV , j eZ . (2.4) 
dxj 

• The formal elliptic Cauchy kernel fl £ A which obeys the following defining equa- 
tions: 

fi-A+*A-ft = O- f^- 1 O P 1 , (2.5a) 
O P - l P Q = OA-'AO. (2.5b) 

in which O is the projection matrix on the central element, i.e. (O • A)j = 5ifiAQj 
and (A ■ 0) iyj = 5 jA it0 . 
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2.2 Symbols of the operators 

The infinite matrix structure described above is most easily explained in terms of the 
symbols of the operators involved. These arise when we consider the action of these 
operators on a basis of monomials, i.e. on an infinite vector c(k) with entries (c(n))j = k? . 
If A denotes a raising operator acting naturally on these vectors by A • c(k) = kc(k) and 
if we allow for the formal expressions P = p(A) — el , 2A • P = p'(A) , with p(/c) = 
p(n\2uj, 2o/) being the Weierstrass p-function with periods 2u,2u/, and where e = p(u>) . 
In this way we obtain a realisation of the operators A and P acting on the vector c(k) as 
follows: 

A • c(k) = fcc(/c) , P ■ c(k) = Kc(k) , 
thus leading to a realisation of the left- index-raising operators: 

(where «i and K2 can be thought of as the symbols for the matrices A resp. *A corre- 
sponding to the rational case). Similarly we can realise their transposeds by their action 
on a corresponding transposed vector *c(k), leading to 

*a ~ C(«2 + ^)-CM-CH = ^ f ( r 2) , ^ - p(« 2 )-e. 

2 p(rc 2 ) - e 

In the above c, C are the well-known Weierstrass functions, and where lu, cj' are the 
corresponding half-periods. Setting to" = —uj — co', then e = p(uj), e' = p(u>'), e" = p(uj") 
are the corresponding branch points of the Weierstrass curve: 

/j. 2 = 4(A — e)(A — e')(A - e") . 

For our later purpose we prefer to write the curve in the following form: 

k 2 = K + 3e + {, (2.6) 
K 

in terms of cordinates (k, K) paramtrised by K = p{n) — e , 2kK = p'{n) , and where 
/ = (e — e')(e — e") . In this way we arrive at the symbols of these operators in terms of 
which we have in particular the realisation of the Caucy kernel: 

o of ^ kl ~ k2 1 - f/( K ^) 

S2 <-» SI{k 1 ,k 2 ) = — — = — — (2.7) 

A 1 - A 2 fcl + «2 

with the identifications: 

h = \ , K i = p(K l )-e (i = l,2). (2.8) 

In order to etsablish some of the results for the Cauchy kernel we need some of the well- 
known addition formulae for the Weierstrass functions. The main identities that we need 
in the context of the present paper are the following ones: 

(C(k + C(«) - CM) 2 = p(« + w) + p(k) + e , 
(p(« + u) - e) (p( K ) - e) = (e - e')(e - e") = / . 
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2.3 Basic relations 

The main object of interest, i.e. the one for which we obtain nonlinear equations, is the 
following infinite matrix U £ A: 

u = c -{l + n-cy 1 , (2.9) 

the components of which are denoted by Uij, € Z). 

In order to obtain the relevant Lax pairs, we introduce also the infinite-component 
vectors u(k) 

u(k) = (1 - U • 0) • c(K)p(k) , (2.10) 

with components (u(K))j, (j E Z), in which p(k) is a scalar function of discrete variables 
n,m E Z as follows 

P(k) = p n , m (k) = (^f) n (P^YpoAk) , (2.H) 



kj \b-k. 

Let us now present the basic equations resulting from this scheme. For the infinite 
matrix U we can derive the discrete matrix Riccati type of relations 

U (a - *A) = (o+A) U-U (O- f l p- 1 OP 1 ) ■ U , (2.12a) 
U-(b - *A) = (6 + A) U- U (O - f l p- 1 O ■ P- 1 ) U , (2.12b) 

together with the algebraic relations 

U *P = P U - U (O A - A O) U ^ (2.13a) 
U ■ *P _1 = P 1 - U + U ■ *P _1 • (O ■ A - *A • 6) ■ P 1 ■ U , (2.13b) 

the latter two equations (which do not involve lattice shifts) being needed for the sake of 
consistency. 

In addition for the continuous case the set of continuous equations can be derived if 
one includes a factor exp(2^^- odd k^Xj) on the r.h.s. of Q2.ll ), yielding: 



_d_ 

dxj 

in which 



U = A j ■ U + U- l A j - U (Oj -f^ 1 Oj -P- 1 ) U , (j odd) , (2.14) 



j'-i 

Oj = ^2(- t A) i -O-A^ 1 ^ . 

i=0 

These equations can be shown to compatible with the discrete relations ( 2.12j ) and lead 
by themselves to a hierarchy of commuting continuous flows. 

In the next section we will derive from the infinite matrix relations ( 2.12j ) and ( p. 14 ), 
together with ( 2,13| ), closed-form partial difference equations (in terms of the lattice shifts 
~ and ^ ) or partial differential equations (in terms of the derivatives with respect to a and 
b) or mixed differential-difference equations, but in terms of specific well-chosen entries of 
the infinite matrix U. Having obtained these equations by means of the infinite matrix 
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structure, subsequently the latter formal tool can be set aside, and the integrability can 
be verified on its own merit. However, in order to obtain additional structures, such as 
for instance the Lax pairs, the infinite matrix structure allows to construct these as well 
in a systematic way. 

Thus, from ( |2.10| ) and the form of p(k), by using the equations ( 2.12| ) and ( 2.1 3| ) together 
with the properties fl2.5a|) of the infinite matrix kernel Q. one can derive the basic set of 
equations for the vector u(k). These constitute the linear discrete relations 



(o-A;)S(k) = (a + A)-U-(0-f t p- 1 -0-P- 1 ) 
(b-k)u( K ) = Ub + A)-U-(0-f t p- 1 -0-P- 1 ) 
as well as the linear "algebraic" relations 

Ku(k) = [P-U (O- A- *A-0)] •«(«) ^ 



u(k) . 



(2.15a) 
(2.15b) 



(2.16a) 

u(«) = [P^ 1 + U ■ t P~ 1 ■ (O ■ A - *A • O) ■ P' 1 ] ■ u(k) , (2.16b) 

P(K) . Similarly, we have 



1 



where K and k are related through the elliptic curve k 2 
linear differential relations for u(k) as a consequence of ( p,14[ ), namely 

d 



u(k) =A j •u(/c) + fc J u(/«)-t7-(Oj - f l P- 1 Oj ■ p-^-ufa) , (j odd) . (2.17) 



From these equations one can eventually derive all the relevant Lax pairs, both for the 
continuum equations as well as on the lattice. 



Remark: The KdV class of systems by the additional requirement that the matrix C 
and hence U is symmetric under transposition, i.e. 



t C = C 



t U = U 



which is equivalent to saying that its entries obey: Uij 
next section in the selection of the components that will enter the main equations 



Uj^. This will be used in the 



3 Elliptic Lattice System 



Having obtained the basic equations ( [2.12 ) and (2.14), together with (2.13), in terms of 
the infinite matrix U we now derive closed-form equations in terms of its entries. It turns 
out that the following choice of entries leads to closed-form equations by combining the 
relations associated with two different lattice shifts, namely 



u 



Uo,o , 
1 - £/-i,o 



U-2,0 , h = 

W=l + [/oi 



U- 



2,-2 



recalling that Uij = Uji. Using the identification of the entries ( |2.1| ) and the relation ( |2.lD 
between the index-raising operators A and P (and similar relation for their transposed), 
we can derive the following set of relations: 
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a + fh 

U-1,-2 + C-1,-2 

a(u — v) 
a(u; — to) 
*7o,i + U ,i 



[a + it)s + v — w 



1 — vw 



a(h — h) — fhh + ss 

fvh + /s (y_ 2 ,-i + U-2-i) + ^-1,1 + 3es + tm 
/uft + fs(U-2-i + ^-2,-i) + 17-1,1 + 3es + 
/ /iu> — s (Uo t i + ?7o,i) + f/-i,i + 3es + urn 
fhw — s (lJo,i + £7),i) + f/-l,i + 3es + wm 
a(tt — -u) + uu — fss 



and similar relations with a replaced by b and ~ replaced by^. 

Combining the various relations we can extract from them the closed-form 
partial difference equations (|1.2| ) which involves only the variables u, s and w. 
tively, we can derive also a lattice system in terms of h, s and v, which reads 



(3.1b) 

(3.1c) 
(3.1d) 

(3.1e) 

(3.1f) 
(3-lg) 



system of 
Alterna- 



a + b + fh-fh) [a-b + fh-fh 



b 2 + f(s-s) 



s — s^j (v — v) = [(a — fh)s — (b — fh)s\ s 
) (%-v\ = ( a + fh)s + {b- flifs 



(3.2a) 



s (3.2b) 



[s — s 



a- fh+Vj (a + A - ~) = a 2 - P(ss) 
b-fh + ^j (b + fh-^j=b 2 -P(ss) , 



a + fh)s- {b + fh)s 
(a- fh)s+ (b + fh)s s (3.2c) 

(3.2d) 

(3.2e) 



and which is obviously equivalent to the lattice system ( |1.2| ). 

If the parameter / = the curve is degenerate and the first equation ( |1.2a| ) will decou- 
ple, yielding the lattice potential KdV equation ( |Ll|) w hich was at the centre of previous 
derivations, cf. e.g. |I7|-||I^|. Thus, the system (|1.2|) can be viewed as a 2-parameter 
deformation of the latter equation. In what follows we demonstrate the integrability of 
the elliptic lattice system (|L~2j), and discuss that the systems is well-posed from the point 
of view of initial value problems. Obviously, we expect that all known results on finite- 
dimensional and similarity reductions that were performed for the lattice potential KdV 
equation can be extended in a natural way to the elliptic lattice system. 

3.1 Lax pairs 

From the relations ( |2.15| ) together with ( |2.16| ) one can derive now a Lax pair for the lattice 
system (|1.2j). In fact, by taking 



(3.3) 
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the following discrete linear systems is obtained: 

[a-k)4> = L(K)(f> (3.4a) 
(b-k)<j> = M(K)cf> (3.4b) 

with the points (k, K) on the elliptic curve representing the spectral parameter. In ( |3.4 
the matrices L and M are given by: 

a — u + -ffJw 1 — j^ss 

L(K) = I K + 3e - a 2 + fss f ~ \ (3-5) 

1 , _ w , f ~ a + u — 4?ws ' 
+{a — u)(a + u) + j^ww " 



and for M a similar expression obtained from ( |3.5| ) by replacing a by b and ^ by ~. It 
is straightforward to show that the discrete Lax equation arising from the compatibility 
condition of the linear system ( |3.4[) 

LM = ML 

gives rise to the set of equations ( |1.2|) . We observe that the matrices L and M depend 
rationally on K only, and thus we do not seem to have essentially a spectral variable on 
the torus. Nevertheless, the solutions seem to depend essentially on the elliptic curve as 
we shall show below by means of soliton type solutions^. 

The Lax pair form the dual system ([^) is derived in a similar way taking the vector 

for which we have a similar linear system: 

(a-k)^ = JL(K)4> (3.7a) 
(b-k)i) = M{K)iP (3.7b) 

where 



( # u i t/ — 4? + 3e- a 2 + fss 

a—fh + Kvs ~, , _ 

+ (a + fh)(a- fh)+Kvv 

y 1 — Kss a + fh — Ksv J 



(3.8) 



and again for M a similar expression obtained from (3^) by making the obvious replace- 



ments. The two Lax pairs (|3.4|) and (pT7h are gauge-related via the gauge transformation: 



</> = KGij, , G = ( * J_! ^ (3.9) 

In fact, the two systems ( |1.2[ ) and (j^^) are essentially the same system, albeit in terms of 
different variables. We note that the gauge conditions on the Lax matrices, e.g. 

GL = LG , (3.10) 

leads to the relations in the (|3.1|), 

lr This seems to be in line with a recent observation by Bordag and Yanovski in where it was shown 
that even for the Landau-Lifschitz equations one can have a Lax pair with polynomial dependence on the 
spectral variable. 
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3.2 Soliton type solutions 

It is relatively straightforward from the infinite matrix structure to construct soliton so- 
lutions, using the symbol representation presented in the previous section. 
Introducing the N x N matrix M with entries 



Mi i = — rj-r r i > (t,i = i,...,JV) (3-11) 



1 - f/(KiKj) 

where the parameters of the solution (ki,Ki) are points on the elliptic curve: 
k 2 = K + 3e + i- , 

and the vector r = (rj)j=i i ... ) jv with components 

'a + fcj\ n 6 + fci ^ m m0 
6 - h 



(3.12) 



where the coefficients r® are independent of n, m. In this case we can take for the the 
infinite matrix C a finite-rank matrix of the form: 

JV 

C = Y,nc^c Ki , (3.13) 

i=l 

and this leads to the following explicit formulae for the quantities of interest: 

u = e-(l + M) -1 -r (3.14a) 

s = e ■ K" 1 ■ (1 + M) _1 ■ r (3.14b) 

to = 1 + e-K' 1 ■ (1 + M)" 1 -fe-r (3.14c) 

u = 1-e-K" 1 -fc- (1 + M)" 1 -r (3.14d) 

/t = e • if -1 • (1 + M)~ l ■ K~ l ■ r (3.14e) 

in which have employed the vector e = (l,l,...,l) and the diagonal matrices 



K = diag(Ki,K 2 , . . . ,K N ) , k = diag(&i, k 2 , . . . , k 



A' I 



We note that although the dynamics itself (encoded in the wave factors r^) does not 
involve the elliptic curve, the soliton solutions essentially depend on the variables on the 
curve. In fact, it is easily verified by direct calculation that the formulae ( |3.14 ) provide a 



solution to the lattice system Ql.2| ) if and only if the elliptic curve relation holds between 
the parameters fc, and the parameters K{. 

4 Initial value problems on the lattice 

We now address the question of how to define a well-posed initial value problem (IVP) for 
the lattice systems (|1.2[ ) and, equivalently, for the dual lattice system (|3.2|) . 
Rewriting the lattice system ( |1.2| ) as follows: 

a + b + u — u ) (a — b + u — u) = a 2 — b 2 + / (s — s) (s — s) (4.1a) 
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w \ „ / w 
a + u — — ] s — [ b + u 



a — u H s+ o + ti — — s 



a + u — — I — u H j = a — P(ss) 



a — u + — ) s — 1 u + — Is 



^ to \ ~ /, ^ it; . 
a + u — \ s + y> — u-\ ) s- 



(4.1b) 

i s 

(4.1c) 
(4.1d) 

(4.1e) 



where as before P(x) = 1/x + 3e + fx , we have to assign values to the dependent 
variables u, s and w on "initial value configurations" of the lattice points, which through 
the eqs. ( 4.1aD -( [4.1ej ) are updated in an unambiguous and consistent way when moving 
through the lattice. As an inspiration for the lattice case it is simpler to investigate first 
what happens in a semi-discrete limit of the lattice system in order to see clearer for what 
variables initial data have to be provided such that we have a welldefined and consistent 
evolution of the the data. 



Through the continuum limit 5 — > , 5 = a — b, employing the Taylor expansions: 

u — > u + 5u + . . . , s — > J+5's + . . . , w — »• w + Sw + . . . , 

in which the dot ' denotes the derivative w.r.t. a continuous time-variable t, we obtain 
the differential-difference system: 

{2a + u-u){l + u) + f{s- s)s = 2a (4.2a) 
[{a + u)s + (a - u)s]s = (s - s)(w + s) (4.2b) 
[w — w — [a — u)s — (a + u)s]s = (s — s)s(l + u) (4.2c) 

a + u-^j (a-u+^j = a 2 - f-4 + 3e + /ssj (4.2d) 



w ws\ ( w\ / w 



1 + j (a - u + — \ + [a + u Ul + ii) = 2a - — + fss (4.2e) 

S S / \ S / \ S ' s s 



s 



It is easy to check that eq. ( |4.2e] ) is redundant as the substition of u, s, w from ( [4.2a ), 



(4.2b) respectively ( [4.2cj ) leads to a triviality. Furthermore, taking the t-derivative of 



( 4. 2d ) and back-substituting the t-derivatives from ( [4.2a|) - ([4.2d ), making use also of ( 4. 2d ) 



as it stands, we get again a triviality. Thus, the DA-system Q4.2p represents a consistent 
time-evolutionary system, for which a well-posed initial value problem can be formulated 
by assigning at t = all values of of u and s along the one-dimensional chain, and providing 
one value of w = wq, say at n = at t = 0. Essentially this means that we are dealing 
with a 2-component DA system in terms of u and s, for which the value wo acts as a 
background value. This is not really surprising, since we are dealing with a two-parameter 
extension of the potential lattice KdV equation and not of the KdV itself. We expect that 
the dependence of the IVP on the background value wq will disappear if one formulates 
the IVP for the properly defined extension of the KdV lattice itself. 
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Let us now return to the fully discrete lattice system, where we can recognise a similar 
situation. Motivated by the KdV case, cf. [21], it is natural to investigate a local iteration 
scheme (i.e. one for which each iteration step at a given vertex involves initial data at 
only a finite configuration of points) is given on "staircases", as in Figure 2, assigning 
values for u and s on the vertices of this staircase, and considering the discrete-time shift 
to be the map (uj,Sj) i— > (Sj,si) . From the consideration of the semidiscrete case we 
expect that we need in addition only one "background" value wq at a specific point on 
the staircase. Indeed, setting the IVP up in this way it can be shown by straightforward 
computation that it is well-posed. In fact, from the three first equations in ( |4.1|) , namely 
( 4.1a| ), Q4,lb| ) and ( |4. lc| ) , one can solve u, s, and w in a unique way, given the values of 
the other variables u, u, u as well as s, s, s and w, w, w. The equations ( 4. Id ) and ( 4.1e ) 
link the variables w and w to w and to the u,s- variables at the relevant lattice points. 

Thus, it remains to be verified that the shifted forms (4. Id) and (4.1e) trivialise through 
back-substitution of w which was already obtained, and this is readily done. Furthermore, 
it is easily checked that the two ways of calculating w from either ( 4,ld| ) followed by 

(4,le) , or from (fl.lej) followed by (4.1 dj) are consistent. 



w 



u ,s 



uo,s |_ 




U3, S 3 



Figure 2: Assignment of initial values on staircases in the lattice. 



We may conclude from these considerations that the coupled system (|4.l| ) is in effect a 
system for for u and s, with one of the eqs. (e.g. (f4.1c ) being redundant, whilst the 
eqs. for w (eqs. ( 4.1d| ) and ( 4.1e[ )) are consistent with the system ( |4.1a| )+( 4.1b| ) involving 
effectively only a background value (or "forcing") through the variable wq. 
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5 Associated Continuous Systems 

As was demonstrated in the past for the lattice systems studied in |L7], [L4|], there 
exist many compatible continuous systems associated with them. These form, in fact, the 
continuous symmetries for the lattice systems (whilst the lattice systems constitute the 
discrete symmetries for the corresponding continuous flows). We will give here a few of 
the simplest of such associated continuous flows for the purpose of identification of the 
associated lattice system. 

Starting from the continuous relations ( [2.14 ) we can derive now a set of partial difer- 
ential equations which are compatible with the discrete equations. Taking into account 
that the derivatives w.r.t. the even time-flows x 2 j (J € Z) all vanish by construction (as 
a consequence of the infinite matrix U being symmetric), we concentrate on the first two 
nontrivial time-flows in terms of x := x\ and t := x%. Thus, from the basic relations ( [2.14 ) 
for j = 1, 3, which can be cast into the forms: 

U x = A-U + U- *A-J7- {O - f l p- 1 -O ■ P- 1 ) U (5.1a) 

U t = \u xxx + hl x - (O - f^- 1 O ■ P- 1 ) -U x (5.1b) 

together with the relation ( |2.13j ) we And the following set of relations in terms of the 
objects defined in section 3: 

2U ,i = u x + u 2 -fs 2 , fh-u= Sx + V ~ W (5.2) 

s 

o o 1 — VW 

217-1,-2 = h x + s 2 - fh 2 , [/!,_! = 5.3 

s 

1 Till} 

-v x = + 3es + v(u + fh) + fs(h x + s 2 - fh 2 ) (5.4) 

s 

Wx = — + 3es + w(u + fh) — s(u x + u 2 — fs 2 ) (5-5) 

s 

for the ^-derivatives, and after eliminating higher powers of A and *A in the infinite-matrix 
relations for the t-derivatives we obtain: 

+ \ u \ ~ \U\ (5.6a) 
3 3 

f -jS x u x - -fh x s x (5.6b) 

+ |*2 - ^ hi (5.6c) 
3 3 

f -v x u x + -fs x (U-l- 2 )x 

3 3 

+ ^s x (U 0> i) x - -fw x h x . (5.6e) 

In order to get the lowest-order system, we eliminate the dependent variables h and v 
and obtain a system solely in terms of u, s and w, which can be simplified further by 
introducing the quantity A = —u + w/s. Thus, we obtain finally the coupled systems of 
nonlinear evolution equations: 





1 


Ut = 


~^v xxx 




1 


St = 


~^xxx 






ht = 


h 

^"'XXX 




1 


v t = 


~7V X xx 




1 


w t = 


~^V) XXX 



(5.6d) 



1 3 

St — ~~^Sxxx 2^ 1 



s z s 2 s 



(5.7a) 
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A t = ± Axxx -^A 2 A x + ^A x (± + 3e + fs^+^(± + 3e + fs^ .(5.7b) 

Although by the introduction of the variable A we have effectively eliminated u and w, 
we mention the following relation among the original variables: 



which can be thought of as the continuous analogue of ( |l.2d ). The coupled system 



is a continuous analogue of the lattice system Ql.2| ), which can be obtained from it by a 
(rather subtle) continuum limit. This system is integrable in its own right in the sense of 
admitting a continuous Lax pair of the form: 



k — u + -j^sw 1 — -J^s 2 



K + 3e - u x - u 2 . f 6 (5.9a) 

2 f 2 k + U- jrWS 



+ 



(U0,l)x Ux 
(Ul,l)x -(U0,l)x 
f ( (1 - Vw)^f +V X W VS x -SVa 



K y (i-vw)y± -w{^) x (vw-i)s*-v x w (5 ' 9b) 

with substitutions: 

uo,i = ^{u x + u 2 - fs 2 ) , (ux t i) x = -(^0,1)^ + u(u 0A ) x - fsw x . 

We have not been able to trace the system of equations (|5.7| ) in the literature, but we 

1 is. 

2 B 



note the amusing fact that if we set A = i — , then the system reduces to 



' ^ s x s xx 3 s x 3 ( \ ,. 2 



st = ^xxx ~ ~^f- + -gf + + 3e + f**j (5.10) 

which is the famous second modified KdV (M 2 KdV) equation, first found by Calogero 
and Degasperis, cf. 0, and also E3, B. However, we must point out that the coupled 
system (|5.7| ) is richer than the M 2 KdV equation, which can be verified immediately by 
reconsidering the soliton solutions ( |3.14| ) presented in section 3.2. In fact, taking the 
factors r% of ( |3.12| ) to include factors of the form exp(2kiX + 2kft) , it is easily verified 
that the formulae ( |3.14 ) provide an explicit solution of ( |5.7| ) with A = —u + w/s. However, 



it is readily checked that this does not provide a solution of ( 5.10| ), nor does the relation 
A = s x /(2s) hold for these soliton solutions! 



It is known that the third modified KdV (M 3 KdV) equation, cf. [11] and also [22], is 



expressible in terms of elliptic functions, cf. also p4J[. However, that equation seems to 



have quite different structure from (|5.7[) , and its soliton solutions, [11|, do not seem to 
involve an elliptic curve. Thus, we are tempted to believe that the coupled system fl5,7| ) 
and its lattice counterpart (^^) forms a genuine 2-parameter "elliptic" extension of the 
continuous and lattice potential KdV equations. 



On a two- parameter extension of the lattice KdV system 



15 



6 Discussion 

In this paper we have presented a general scheme to obtain integrable systems (i.e. par- 
tial difference equations) associated with elliptic curves. These systems constitute a 2- 
parameter deformation of the lattice systems that were investigated in the past. For the 
case / = 0, i.e. when the elliptic curve generates into a rational curve, we recover the usual 
lattice KdV system that was studied from various perspectives during the last decades. 
One relevant recent discovery in which these systems played an instrumental role, is that 
the KdV hierarchy -when formulated appropriately- is rich enough to contain the full 
Painleve VI equation as similarity reduction. This slightly surprising observation was 
made in a recent paper |l8| and subsequently it was shown in [16] that the entire KdV 



hierarchy can be encoded in one single nonautonomous PDE which is also reducible to the 
full PVI equation with arbitrary parameters. This demonstrates that these systems are 
in some sense universal. The programme that now unfolds on the basis of the results of 
the present paper, is to develop similar systems associated with an elliptic curve, which 
would then constitute 2-parameter deformations of the systems just mentioned. In that 
way one expects that these systems under appropriate similarity reductions reduce to an 
"ellipticised" version of PVI, possibly the equations that were derived in the 1970's by K. 
Okamoto from isomonodromic deformation problems on the torus, cf. [p0| . 

We should mention that there exists also an alternative way to extend the lattice 
systems of KdV type such that there is an underlying elliptic curve. In fact, V. Adler 
discovered in a lattice version of the famous Krichever-Novikov equation, cf. |jlC[ |. The 
resulting lattice system has the following form: 

kQUUuu — k\ (^uuu + uuu + uuu + unuj + k<i (uu + u-uj 

— &3 (uu + uu ) — &4 (uu + uu) + k$ ( u + u + u + u ) + k$ = , (6-1) 



where the coefficients ko, ■ ■ ■ , ke are parametrised through elliptic functions. The main dif- 
ference between ( |6.1[ ) and ( |1.2[ ) resides in the fact that the lattice parameters for the Adler 
equation are themselves points of the elliptic curve, namely with coordinates (p(a), p'(a)), 
(p(/3), p'(/3)), whereas this is not the case for the system (p..2|). Furthermore, the Lax pair 



for (6J), which was constructed in [13], has the spectral parameter living also on the ellip- 
tic curve, whereas this is not evident for the Lax pair for (|L2j). Nonetheless, the explicit 
formulae for soliton solutions which were discussed in subsection 3.2, demonstrate very 
clearly that they contain parameters (ki,Ki) which live on the elliptic curve. Thus, this 



seems clear indication that the elliptic curve ( |2,6D is essential, and that the lattice system 
( |1.2| ) constitues a genuine 2-parameter extension of the lattice (potential) KdV equation 
( |1 . 1| ) . A further comparison of between the two lattice systems is still needed as well as a 
clear identification of the underlying continuous system ( |5.7| ) . 
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